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A program was initiated to improve the flow quality 
through the following steps: 

1) CFD modelling of the existing flow, 

2) Development of possible nozzle contours, 

3) Evaluation of these nozzle contours 

4) Remachining of the nozzle, and 

5) Validation testing. 

Mitsubishi Heavy Industries (Mi ll), Aero Systems 
Engineering (ASE), and NASA Langley participated in 
this effort with NAL. MHI did the contour 
measurements and remachining; MHI and NASA LaRC 
used their CFD capability to model the original test 
results and evaluate the new contours; and ASE 
collaborated to develop the improved nozzle contours. 

ANALYSIS - EXISTING CONTOUR 

The initial test data shown in Figure 1 was expanded to 
include data at several other axial stations in the test 
section. As shown in Figure 2, the data indicate 
pressure waves traversing the test section and the 
effects of other waves converging on centerline. The 
spread of Pitot pressure ratio is + 7%; this implies a 
Mach Number variation of + 1.4%. This profile was 
shown to be nearly symmetrical. 

Through a preliminary CFD analysis these waves were 
speculated to originate from the nozzle wall as shown 
in Figure 3 

Figure 1 Pitot Pressure Profile - Center of the Test 
Section 
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SUMMARY 

As a result of CFD analysis and remachining of the 
nozzle, the flow quality of the Mach 10 Hypersonic 
Wind Tunnel at NAL Chofu, Japan was improved. 
The subsequent test results validated the CFD 
analytical predictions by NASA and MHI. 

INTRODUCTION 

The 1 .27 meter diameter Mach 10 nozzle at the 
National Aerospace Laboratory (NAL) at Chofu, 
Japan was designed in the early 1990s for a nominal 
operating condition of 1070 K and 6 MPa. This 
nozzle was designed and manufactured in accordance 
with industry standards at the time. The flow quality 
(see Figure 1) of this project was acceptable at the 
time l ; but since then, the standards of the industry 
have been tightened 2 , and the Chofu design was not 
up to recent w'orld-class performance criteria 3 . 
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Figure 2 Pitot Pressure Profiles - Various axial 
stations 
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Figure 3 Nozzle Sections and Possible Wave 
Pattern 


While one wave could come from the joint between 
Expansion 1 and Expansion 2, the others were not so 
obvious. Those others could be reflections of waves 
that originated further upstream or they could be 
waves generated on the nozzle wall at the location 
shown. Therefore, the wall contour was precisely 
measured at many stations so small disturbances 
could be detected. The results of these measurements 
showed that the contour was, in general, a shape 
typical of hypersonic nozzles but the slope and 
curvature were not fully smooth as shown in Figure 4 
for the throat and initial expansion. 



The complete measured contour data was placed in 
both the NASA-Langley 4A6 and MHI w CFD 
models and the predictions are shown in Figure 5 and 
Figure 6 respectively. In each case, the general 
character of the non-uniform profiles is predicted. 


(The differences that do exist are not deemed to be 
significant because a prediction of such steep gradients 
is subject to exact measurements of all peaks and 
valleys and exact duplication of flow properties over 
the entire flow field. While the coordinates of the 
contour were measured in fine detail, there is the 
possibility that some peaks/valleys were not recorded 
and/or that some gas properties are not fully modeled.) 



Figure 5 Comparisons of Pitot Pressure Profiles 
with CFD - NASA Langley 




Figure 6 Comparisons of Mach Number Profiles 
with CFD - MHI 

OTHER TOOLS FOR ANALYSIS 

Along with the Pitot pressure profiles, the usual contour 
maps of Mach, pressure, velocity, etc. were also 
obtained (see Figure 7 for an example). The origins of 
waves were not easily interpreted from these usual 
contour plots and the parameter "GradP/P” was 
developed. These contours are shown in for the original 
nozzle. As shown in Figure 8, many waves that come 
from the wall can be seen. Some of these waves get 
cancelled through various interactions and a few waves 
enter the test section. Two of these waves generally 
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agree with the two pressure waves that were identified 
in Figure 2. The strong gradients on the centerline are 
also obvious in this Figure 8. Contours could be 
analyzed more easily on the basis of this CFD 
parameter. 



Figure 7 Typical Mach Contour 



Figure 8 Contour of GradP/P for Original Nozzle 


CON FOUR SMOOTHING 
Initial Smoothing 

Based on the successful validation of the CFD 
models, new contours were developed in an attempt to 
obtain more uniform flow. One of the criteria for 
these new contours was that most of the existing 
nozzle had to be re-used in order to minimize cost. 

The sections of the nozzle are shown in Figure 3. 

The first attempt was to smooth just the throat block; 
that is, the initial expansion of the nozzle. This is a 
small section and could be made new' for relatively 


low cost. The results are shown in Figure 9. Note that in 
both the cases that were studied, major non-uniformities 
still remain. In this case, it does not seem possible to 
achieve uniform flow with modification of only the 
throat block. 
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Figure 9 Pitot Pressure Profile for Two Cases of 
Throat Block Smoothing - Compared w ith Original 

Final Smoothing 

Further attempts to obtain uniform flow consider the 
entire nozzle, not just the throat block. The onginal 
contour was based on the Method of Characteristics 
plus Boundary Layer correction (MOC+BL) and was 
then smoothed using cubic spline methods. For the new 
coordinates, a revised smoothing method was used 
based on polynomial equations of nozzle radius as a 
function of nozzle axial location. The nozzle was 
divided into axial sections for the curve fit procedure. 

In each section a 5 th order polynomial equation w'as 
applied. The polynomials were of the 5 th order since 
this allowed the necessary inflections of the curvature. 
(For reference, a cubic equation has a linear second 
derivative, i.e. curvature, and use of several cubic 
equations would have a discontinuous third derivative, 
not the smooth charaterislic needed for a nozzle.) The 
lengths of the sections were adjusted to provide a good 
fit of the existing coordinates. Five sections were 
needed as shown in Figure 10. 

At the match point between each section, the 
coefficients were adjusted until there was a match for 
nozzle Radius, Slope, and Curvature. The section from 
about 2500 mm to 5500 mm is covered by a single 
equation. This is a subtle but important point; any wave 
•originating in this range goes directly to the core of the 
test section, virtually without attenuation. This section 
must be a very smooth and this is best provided by a 
single equation. 
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Figure 10 Sections of Nozzle for Curve Fit 
Procedure 


Further, since the MOC+BL method is step-by-step 
process that uses assumptions and averages to 
compute the downstream points, it is possible that the 
best coordinates might be slightly different than the 
computed values. Therefore, exact matching of the 
MOC+BL coordinates was not a priority in this 
smoothing process. It is more important to have a 
continuous third derivative (smooth curvature) than to 
exactly match the MOC+BL computed points. 

One criteria for this nozzle was that the existing 
surface had to be machined at least 0.5 mm so that the 
new contour could be made with confidence. This 
allowance was needed for the tolerances of machine 
setup, concentricity, and a full removal of the existing 
coating. This was accomplished by rotating the new 
contour about a point upstream of the throat. (The 
rotation was only 0.0005 radians; additional 
smoothing was done in the throat area to maintain the 
throat size.) 


NEW CONTOUR 


The above method was used to generate several 
nozzle contours that might improve the flow 
uniformity. Each contour was evaluated using the 
CFD already validated to the test data as mentioned 
above. 

The final selected change in contour is shown in 
Figure 1 1 . This shows the slight modification that was 
applied to the original nozzle. Note that, exclusive of 
the general downward trend (resulting from the 
rotation mentioned above), the contour change is only 
about +0.25 mm. The peaks in the line at about 3400 
and 5000 mm indicate compression comers (as 
compared to the new contour) and they nearly match 
the origin of waves suggested in Figure 3. At each of 
these two locations, the slope is changed by about 1.6 
mm/ 1000mm (~0. 1 deg total change); this highlights 


the need to limit the size of slope variations, especially 
in the compression mode. Note that the sharp 
compression comer near 700 mm may have contributed 
to the problem on centerline (after several reflections); 
for reference see Figure 8 that shows several strong 
waves originating in this area. 



Figure 11 Change in Contour 


MACHINING RESULTS 


The results of the machining to the new contour are 
shown in Figure 12. The contour was held to within 
nearly +,- 0.025 mm of the design values. The resulting 
compression comers near the joints at 2500 and 5800 
mm are less than 0.01 deg. There is a step, about 0.03 
mm high, near 2500 and a similar gap near 4 1 GO mm. 
As will be shown later, the flow quality is quite good 
despite these deviations. 



Axid Location - mm 


Figure 12 Results of Machining to New Contour 
TEST RESULTS 

Pitot Pressure 

Several improvements were made in the measurement 
apparatus prior to conducting the final tests. The 
number of Pitot probes was increased and the span 
extended so that measurements are at 10 mm transverse 
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spacing and the measurements now extend into the 
boundary layer. The test data was taken at more axial 
stations, with 50mm increments in the axial direction 
covering 750 mm, in both the horizontal and vertical 
planes including tests with off-axis placements (again 
both in the vertical and horizontal orientations). The 
data is summarized in Figure 13; note that the spread 
of Pitot pressure within the core +350 mm around 
centerline has been reduced from +7.5% to +1.5%. 
This implies that the Mach Number variation is now 
been reduced from + 1 .5% to ± 0.3%. 


Summary of Test Data After Modification 



Figure 13 Final Test Results 

Flow Angle 

The flow angularity was measured using a rake with 
conical probes. The measurements cover 0.72 m in 
width with 7 conical probes spaced 0. 12 m. These 
tests were made in 10 axial locations, every 0.1 m 
from nozzle exit (X = 0) to 0.9 m downstream, with 
the rake in both the vertical and horizontal position, 
totaling 130 points. With the rake in a horizontal 
position, the distribution of flow angle was obtained 
in 7 vertical locations (Z = 0, ±0.12, ±0.24, ±0.3 m) at 
3 longitudinal locations (X= 0, 0.25, and 0.5 m), 
totaling 147 points (some duplicating the longitudinal 
positions). 

The flow inclination angles were from -0.02 degrees 
to 0. 1 degrees on the centerline and within ±0.2 
degrees in the whole core flow region. An example of 
the result at the nozzle exit cross section is shown in 
Figure 14. The flow is nominally axisymmctric with 
slight inward flow at this location. 

Centerline Mach Number 

Further evidence of the flow quality achieved in this 
effort is shown in Figure 15 (taken from Reference 
10). 'This shows the comparison of the centerline 
Mach number in this wind tunnel as compared to the 
same parameter in the AEDC Tunnel C. The long 


length of nearly constant Mach number in the NAL 
tunnel is apparent. 


Nozzle exit cross section 
(1.27m Dia.) 



Figure 14 Flow inclination angle distribution at 
nozzle exit cross section 



Figure 15 Centerline Mach Number Variation - 
Compared with AEDC Tunnel C (from Reference 
10 ). 

CONCLUSION 

An international team has studied the flow quality of 
the NAL Mach 10 nozzle. The participants used both 
experimental data and computational methods to 
identify concerns and to modify the nozzle contour. The 
results with the new contour show that the flow quality 
has been significantly improved. As a result the NAL 
Mach 1 0 wind tunnel can now serve the aerodynamic 
community as a world-class testing facility. 
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